Hole-depletion of ladders in Sri4Cu2404i induced by correlation effects 
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The hole distribution in Sri4Cu2404i is studied by low temperature polarization dependent O K 
Near-Edge X-ray Absorption Fine Structure measurements and state of the art electronic structure 
calculations that include core-hole and correlation effects in a mean-field approach. Contrary to all 
previous analysis, based on semi-empirical models, we show that correlations and antiferromagnetic 
ordering favor the strong chain hole-attraction. For the remaining small number of holes accommo- 
dated on ladders, leg-sites are preferred to rung-sites. The small hole affinity of rung-sites explains 
naturally the ID - 2D cross-over in the phase diagram of (La,Y,Sr,Ca)i4Cu2404i. 

PACS numbers: 78. 70. Dm, 71. 15. Mb, 74.78.-w, 74.72. Gh 



I. INTRODUCTION 

The quasi-one-dimensional spin chain and ladder 
(La,Y,Sr,Ca)i4Cu2404i compounds have attracted con- 
siderable interest since the discovery of a quantum criti- 
cal phase transition in their phase diagraniji These com- 
pounds are the first superconducting copper oxide ma- 
terials with a non-square lattice. They are composed of 
alternately stacked chain and ladder planes. The chains 
are made up of Cu02 linear edge-sharing CUO4 squares 
and the ladder planes consist of two zig-zag strings of 
corner-sharing CU2O3 squares (see Fig. [1]). The parent 
compound Sri4Cu2404i is naturally doped with six holes 
per formula unit [fu). 

Determining exactly how the holes are distributed in 
the system is hindered by the complexity of the crys- 
tal structure of Sri4Cu2404i and electron correlation 
effects i^i^ Room temperature optical conductivity;^ Near- 
Edge X-ray Absorption Fine Structure (NEXAFS)^ 
X-ray emission spectroscopy^ and Hall coefficient 
measurements^ on Sri4Cu2404i estimated that at least 
five holes per fu reside in chains and at most one in 
ladders. High density of holes in chains is compatible 
with the low temperature Tc ~ 200 charge ordering 
with the 5-fold chain periodicity accompanying the an- 
tiferromagnetic (AF) spin-dimerization, observed by in- 
elastic neutron scattering 1^"— In ladders, a gapped spin- 
liquid and the charge density wave (GDW) appear in 
the ground state, in the spin and charge sectors, respec- 
tively. Whereas the spin-liquid and its collective spin ex- 
citations (triplons) are well understood theoretically and 
experimentally)^"— the low hole density in the charge 
sector is at variance with the 5-fold periodicity in ladder 
cell units (Wigner hole crystal) observed by low temper- 
ature resonant X-ray diffraction measurements at the O 
K edge if a 4ki? GDW picture is assumedJ^ On the other 
hand, a revisited interpretation of NEXAFS spectrapi^ 
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FIG. 1: (Color online) 3D view along the c crystal axis of four 
316-atom AF unit cells of Sri4Cu2404i. The structure is com- 
posed of alternately stacked chain layers, constituted of CUO2 
linear edge-sharing CUO4 squares, and ladder planes, made 
up of two zig-zag strings of CU2O3 corner-sharing squares. 
They are parallel to the (a,c)-plane and separated by Sr 
atoms. Cu/O/Sr atoms are colored in red/blue/yellow. 



claimed a distribution of 2.8 holes in ladders and 3.2 
holes in chains, which satisfactorly explains the GDW 
in ladders, while it fails to explain AF dimer order in 
chains. 

The apparently contradictory results are mainly re- 
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lated to the lack of a suitable theoretical model. All 
published NEXAFS analyses are based on semi-empirical 
models that do not take explicitly into account the com- 
plex interplay between spins and holes in ladders and 
chains. Thus the basic question of where the holes reside 
in Sri4Cu2404i has remained unresolved up to now. 

The hole distribution between chains and ladders is the 
key to stabilization of interdependent electronic phases 
in these two subsystems. It is a prerequisite to under- 
standing the low-energy physics of the system and such 
open issues as the occurrence of Zhang-Rice singletSfii 
the evolution of spin order and the origin of the su- 
perconductivity in the Ca-doped compoundi^ We bridge 
the gap between theory and experiment by performing 
low temperature electron yield polarization dependent O 
K pre-edge NEXAFS measurements and state-of-the-art 
electronic structure calculations. In the theoretical mod- 
eling we consider the full 316-atom AF unit cell, include 
the core-hole effects, core level shift and correlations in a 
DFT-I-U framework. This method has proven to be fairly 
successful in reproducing X absorption spectra in corre- 
lated metals with well localized orbital o^^i^^ . We show 
that, contrary to previous claims, correlation effects and 
AF order stabilize holes on chains, induce hole-depletion 
in ladders, where rung sites become less populated than 
leg sites. 



II. EXPERIMENTAL AND CALCULATION 
DETAILS 

High-quality single crystals of Sri4Cu2404i were grown 
by the traveling-solvent floating-zone method and char- 
acterized by X-ray diffraction measurements. A well ori- 
ented sample was cleaved in-situ along the (a,c) plane 
under a pressure of 10~^° mbar. The O K edge NEX- 
AFS measurements were performed using the helical un- 
dulator Dragon beam line ID08 at the European Syn- 
chrotron Radiation Facility (ESRF) in the total electron 
yield mode with 200 meV resolution. The incident light 
was normal to the sample surface and its polarization 
was parallel to the a/c sample axis. 

NEXAFS calculations are performed by the 
XSPECTRAi^-^ code based on Density Func- 
tional Theory (DFT)2i and all-electron wave function 
reconstruction)^ Correlation effects are simulated in a 
mean field DFT-t-U framewor k^^'^^ using a Hubbard on- 
site energy of U^rf = 10 eV for copper and Upp = 4 eV 
for oxygen. The chain and ladder atomic positions 
were taken from Gotoh et ali^ without imposing any 
superstructure modulations, either along the chains or 
along the ladders. We performed a specific NEXAFS 
calculation for each of the five inequivalent oxygen sites 
existing in the structure, two in the ladders, and three 
in the chains. We find a negligible core- level-shift. 
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FIG. 2: (Color online) Polarization dependent O K NEXAFS 
pre-edge intensity for E|[a and E||c measured at 150 K.). Sub- 
script e is for experiment. 
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FIG. 3: (Color online) Temperature dependence of the O K 
NEXAFS pre-edge intensity for E||c showing that for this po- 
larization spectra change below and above the charge ordering 
temperature (« 200 K). The inset shows that the E||a spectra 
temperature change is negligible. 



III. RESULTS AND DISCUSSION 

A. NEXAFS measurements 

The low temperature (150 K, see Fig.© O K NEXAFS 
pre-edge spectra has two distinct structures, Ae, centered 
at 527.9 eV and Be, at 529.5 eV (e is for "experiment" ) , in 
agreement with previously published room-temperature 
spectrai^i^ii^ The shape of the two structures is depcn- 
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dent on the polarization of the incident photon: for E||a, 
Ae is stronger than Be, while for E||c (direction along the 
chains and the ladders) it is the opposite. For E||c struc- 
ture Ae has a second, well resolved peak at 528.5 eV, 
A'e. The intensity of this structure is only weakly af- 
fected by the charge ordering transition (see Fig. [3]). The 
E||a spectra are even less affected at Tp (see Fig. [3] in- 
set) . The attribution of the features in these spectra has 
so far been very controversia h^i^^'^^ 



B. Calculation of the NEXAFS spectra 

Theoretical modeling using the AF unit cell with 5- 
fold periodicity in chains is shown in Fig. |4l The two, 
Ae and Be, structures are well reproduced in the calcu- 
lations (labeled Ac and Be where c is for the calculated 
spectra), although their separation is only 0.7 eV, com- 
pared to 1.6 eV for the experiment. This discrepancy 
is attributed to the well-known underestimation of the 
Hubbard gap in DFT-I-U simulations. The agreement 
between theory and experiment is particularly good for 
the E||c polarization as in the c-direction the charges arc 
less localized. For E||c (i) the Be cross section is stronger 
than the Ac cross section, and (ii) the high-energy shoul- 
der A'c is clearly visible. For E||a polarization, the Be 
structure has smaller cross section than the Ac peak. 

The two insets in Fig. H] show the calculated cross sec- 
tion contribution of chain (C), ladder-leg (L) and ladder- 
rung (R) oxygens to the NEXAFS spectra. The con- 
tribution from ladder-rungs is small or negligible for all 
features except for the Be peak in the E||a geometry. 
This result is in line with the analysis by Niicker et al<^ 
and is in strong disagreement with the results of Rusydi 
et al.— where the contribution of ladder rungs was over- 
estimated. Our results demonstrate that ladder-rungs 
are hole-depleted in Sri4Cu2404i. 

For both polarizations the low energy Ac peak is 
equally composed of chain and ladder-leg states, with 
a minor contribution from ladder-rungs in the E||a ge- 
ometry. Its high energy A'c shoulder is composed of 
both ladder-leg and chain contributions. This result is 
in stark disagreement with all previous claims based on 
semi-empirical NEXAFS analysis^ii^ that attribute the 
A'c structure solely to ladders. It points to the fact that 
the O K edge resonant diffraction results, performed at 
the photon energy corresponding to A'c, could be related 
to the modulation reflections due to the chain-ladder lat- 
tice mismatcb2£. The Wigner-hole crystallization analy- 
sis of the O K edge resonant soft X-ray scattering data 
in Ref. [T^ now becomes questionable. 

Calculations using the AF unit cell with 4-fold period- 
icity in chains^! is shown in Fig. [5j The two insets show 
the calculated cross section contribution of chain (C), 
ladder-leg (L) and ladder-rung (R) oxygens, equivalent 
to these in Fig. The experimental results are less well 
reproduced compared to the 5-fold chain-periodicity unit 
cell. First, the Ac - Be energy separation is smaller. Fur- 




I 

-10 12 3 

Energy (eV) 



FIG. 4: (Color online) Calculated O K NEXAFS pre-edge 
cross section on the AF unit cell with 5 fold periodicity in 
chains, for E||a and E||c. Insets : contribution of chain (C), 
ladder-leg (L) and ladder-rung (R) oxygens for E||a (upper 
inset) and E||c (lower inset). Zero energy corresponds to a 
photon energy of 527.4 eV. Subscript c stands for calculated. 
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FIG. 5: (Color online) Calculated O K NEXAFS pre-edge 
cross section on the AF unit cell with 4 fold periodicity in 
chains, for E||a and E||c. Insets : contribution of chain (C), 
ladder-leg (L) and ladder-rung (R) oxygens for E||a (upper 
inset) and E||c (lower inset). Zero energy corresponds to a 
photon energy of 527.4 eV. 



ther, for the E||c geometry Ac and Be are almost equal 
in intensity, contrary to the experiment. Finally, A'c ap- 
pears as a low energy shoulder to structure Be and not 
as a high energy shoulder to Ac. For these reasons, the 
4-fold chain-periodicity unit cell can be eliminated as a 
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model for chain spin ordering. 
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FIG. 6: (Color online) (upper panel) Total DOS of 
Sri4Cu2404i without a core-hole, calculated on the AF unit 
cell with 5 fold chain periodicity. The Fermi level (Ef) corre- 
sponds to eV. (lower panels) Local chain- and ladder-DOS 
projected to Cu (red line) and O (blue lines). The second 
dotted line indicates where Ef would be if the system had 
6 additional electrons per fu, i.e. no holes (Ef' = Ef+AE, 
where AE = 0.31 eV corresponds to additional 6 x 4 = 24 
electrons) . 



C. Density of states 

To understand the role of core-hole effects we also per- 
formed density of states (DOS) calculations for the sys- 
tem without a core- hole (see Fig. 0]). Contrary to previ- 
ous theoretical works r^"— performed on a minimal cell 
without AF ordering, we find the occurrence of a 0.25 eV 
gap appearing in ladder states. This finding is in agree- 
ment with the insulating behavior of the ladders^"— and 
points out the need of including the full AF crystal struc- 
ture in the simulation. The occurrence of a 0.25 eV Hub- 
bard gap in ladder states results in a strong hole deple- 
tion of this sub-system. Indeed, we find that only about 
0.4 holes reside in ladders, while the chains accommo- 
date about 5.3 holes per fu. A more detailed analysis 
of hole distribution is presented in Table HI The mean 
value of the probability of accommodating a hole on one 
atom is schematically presented in Fig.|8]by the thickness 
of the contour surrounding atoms of each group (chain, 
ladder-leg, ladder-rung). 

The role of correlations is even more evident when an- 
alyzing hole and spin order along the chains in the para- 
magnetic DFT calculations and in the DFT including U 
and AF ordering. When correlation effects are switched- 
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FIG. 7: (Color online) Chain site dependent hole (spin) den- 
sity : (a) Chains have natural 5-fold hole periodicity in the 
paramagnetic cell without correlations ; (b) When the 4-fold 
chain spin periodicity is imposed, the 5-foId hole periodicity 
is destroyed ; (c) In the case of the 5-fold spin periodicity, the 
natural hole periodicity is preserved, on oxygen and as well 
on copper sites. 
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TABLE I: Calculated fraction (n) of 6 holes per fu residing in 
chain-, ladder- and in Sr-inter-layers of Sri4Cu2404i; number 
of atoms in each group {Ng) ; distribution of holes between 
Cu- and 0-sites in chains and ladders (rig). 



on the number of holes in the ladder is decreased by a 
factor of five. Hole- depletion of ladders is thus a pure 
correlation effect. 

The distribution of more than 5 holes in chains and 
less than 1 in ladders is compatible not only with the AF 
unit cell with 5-fold chain periodicity, but also with the 
AF model with 4-fold chain periodicityi^ But the latter 
model turns out to be irrelevant since in this case the 
NEXAFS calculations could not reproduce the experi- 
mental data (see Sect llll B|) . Moreover, the 4-fold chain- 
spin periodicity destroys the natural 5-fold hole ordering, 
while the 5-fold spin periodicity preserves it, as shown in 
Fig 13 We would like to stress once more that taking 
into account the proper long-range AF order is crucial 
to achieve a correct description of O K-edge NEXAFS 
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CDW ground state is then suppressed in favor of the 2D 
superconductivity under pressureJ^ 
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FIG. 8: (Color online) The thickness of the contour at 0/Cu 
sites symbolizes the mean hole density per site in chains (up- 
per panel) and ladders (lower panel) as determined by ab- 
inition calculations. 



spectra. 

It is worthwhile noting that our DOS is in good agree- 
ment with experiments as we find the occurrence of a 
strong -2 eV band with mixed chain/ladder contribu- 
tions and large Cu-ladder weight, in accordance with 
Angle Resolved Photoemission Spectroscopy (ARPES) 
measurementsj^ Furthermore the two ladder-DOS struc- 
tures at -0.25 eV and -0.1 eV have also been observed in 
a recent ARPES^ experiment and were identified as the 
quasi-lD underlying Fermi surface. 

Finally, the tendency of hole-depleting ladder-rungs 
found in our work, provides a satisfactory expiation for 
the 2D to ID cross-over when the number of holes in 
the system is reduced, like in under-doped compounds^ 
or at low temperature where the less localized ladder 
sub-system suffers from the migration of holes to chains 
(back-transfer) as reported in the NMR<^ On the fully 
doped side we suggest that the Ca doping has a twofold 
effect: it destroys long-range AF order in chains2i22i2^ 
and reduces hole-depletion in ladders. Increased popu- 
lation of the rung sites reinforces the 2D character of 
the ladders. For sufficient Ca-induced ladder-doping the 



IV. CONCLUSION 

In conclusion, our experimental NEXAFS spectra and 
first-principles theoretical modeling including the com- 
plete 316 atom AF unit cell, core-hole and correlation 
effects demonstrate that holes mainly reside on chains 
and thus ladders are hole-depleted. This finding resolves 
long standing controversial interpetations by several au- 
thors based on semi-empirical models. The analysis by 
Niicker et ali^ is in agreement with our findings based 
on an ab-initio approach. It is clear, however, that their 
arguments were based on a too limited description since 
correlations and long range AF behavior consisting of a 
5-fold chain periodicity must be included to obtain the 
corrrect interpretation of the experimental spectra. 

In ladders, majority of holes reside in leg sites, while 
only a tiny minority populates rung sites. Further exper- 
iments and calculations arc under consideration to verify 
our suggestion that the small affinity of rung sites in the 
parent compound explains the ID - 2D cross-over going 
from the under-doped to the fully doped side in the phase 
diagram of (La,Y,Sr,Ca)i4Cu2404i. 

Our work unambiguously answers the out-standing 
question of where the holes reside in Sri4Cu2404i. This 
is crucial to understanding all phenomena occurring in 
the phase diagram of this family of compounds, ranging 
from charge-ordered antiferromagnetism to superconduc- 
tivity and is a prerequisite to all investigations based on 
correlated models of spin ladder and chain compoundsi^ 
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